Temporal Scaling of Interfaces Propagating in Porous Media 



G\ 
On 

S3 
33 
in 

OO 



> 

o 

m 

o 

ON 

C3 



C3 
O 

o 



13 



Viktor K. Horvathi and H. Eugene Stanley 
Center for Polymer Studies and Department of Physics, 
Boston University, Boston, Massachusetts 02215 
(Physical Review E, pp. 5166-5169, 52 1995 Nov.) 

To better understand the temporal behavior of a roughening meniscus driven by capillary forces 
during the imbibition of a viscous fluid in porous media, we measure the height-height autocorre- 
lation function C(t) using a constant driving force. We find C(t) ~ t 13 , with (3 — 0.56 ± 0.03, and 
provide the first experimental evidence for driving force independent temporal scaling behavior of a 
propagating wetting front in the presence of quenched noise. We interpret the value of /3 in terms of 
the possibility that the dynamics may be governed by nonvanishing non-linearity due to anisotropic 
depinning. 



The growth of rough surfaces and interfaces under far- 
from-equilibrium conditions is a common phenomenon 
in nature jjj]. Examples include such processes as vapor 
deposition, crystallization, thin film growth by atomic 
beams, settling of granular materials, and fluid flow in 
porous media. The fluctuations of the interface height 
h(x,t) can be characterized either by their standard de- 
viation a(l,t), or by the height-height correlation func- 
tion C(l,t) = [((h(l + x, t + r) - ft(^,r)) 2 ) XiT ] 1 /2 ) where 
h(x,r) = h(x,r) — {h{x,r)) x . The common belief is 
that cr(l,t) and C(l,t) exhibit the same general statis- 
tical properties. 

Particularly useful is the Family- Vicsek dynamic scal- 
ing hypothesis M, 



C(t)=C(0,t) ~ L a F(t/L a /P), 



(1) 



where L is the system size and !F{y) ~ y 13 for y <C 1, 
T(y) ~ const for y ^> 1. Extensive numerical simula- 
tions of different computer models [jlj predict universal 
values for the exponents (3 and a. Theoretical work is 
based on Langevin-type equations, such as 



dh(i,t) 
dt 



v 



vV 2 h + A e// (V/i) 2 +??(x,<). (2) 



The linear [X e ff = 0) an d the non-linear (X e ff 0) cases 
are referred to as the Edwards- Wilkinson (EW) (|] and 
Kardar-Parisi-Zhang (KPZ) 0] equations. Analytical so- 
lution of Eq. (j^) in d! = 1 (where d' is the dimension 
of x) reveals the same spatial exponent (a = 1/2), but 
different temporal exponents 



f3 EW = l/A, KP z = 1/3. 



(3a) 



If we know these exponents, then we can distinguish be- 
tween linear and non-linear dynamics based on p. There- 
fore it is important to make accurate measurements of /?. 

If the randomness in the dynamics is due to the in- 
homogeneity of the media where the moving phase is 
propagating, then the resulting interface exhibits differ- 
ent scaling behavior. This time- independent "quenched" 
randomness can be described by replacing ?7(x, t) by 



n(jL,h(t)) in Eq. (||). We denote the linear and non- 
linear equations as QEW and QKPZ respectively. Nu- 
merical solution of the QEW equation || produces self- 
affine interfaces with a roughness exponent a in the 
range [0.5 — 1.0], that is tunable with the driving force 
vq ]o]- ^0| , ^2| -[r^] . Different computer models also exhibit 
very scattered results for a 0. However, theoretical con- 
siderations suggest H that Pqew — (4 — d')/A. Numeri- 
cal integration of the QKPZ equation close to the pinning 
transition exhibits temporal scaling with f3 = 0.61 ±0.06, 
in good agreement with the result of simple dimensional 
analysis which predicts Pqkpz = (4 — <f')/(4 + d') ||. 
Hence for d' = 1: 



Pqew — 3/4, Pqkpz — 3/5. 



(3b) 



The fact that scaling of surfaces during roughening in the 
presence of quenched noise exhibits driving force inde- 
pendent temporal scaling also underlines the importance 
of the investigation of the exponent fi, which is the main 
goal of this Letter. 
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FIG. 1. Schematic of the experimental setup. The av- 
erage height h of the meniscus in the wetted paper was 
kept at a constant level. A video camera digitized the 
interface in real time and the computer pulled down the 
paper if h became larger than the predefined value H. 

Specifically, we investigate the growing interface during 
imbibition of viscous liquids in filter paper. The experi- 
mental setup (Fig. 1) contains two vertically-positioned 
40cm x 30cm parallel plexiglass plates separated from 
each other by 3 cm and closed with side walls. The bot- 
tom part of the plates extend into a liquid container from 
which water starts to wet the filter paper. We designed 
the size of this liquid container to be large enough to 
keep the free surface of water at a constant level. All 
height measurements are relative to this reference level. 
The upper part of the cell is closed by polyethylene film 
to prevent a large amount of evaporation. The dry pa- 
per passes through two touching cylinders driven by a 
stepping motor, and enters into the cell through a nar- 
row gap on this film. We use other rollers and a proper 
weight to keep the paper stretched. Stretching of wet 
paper sometimes involves changes in the paper structure 
mostly by elongating the paper strip. This was checked 
by double measurements of the speed of the paper strip. 
We measured the speed of the paper both at the driving 
rollers at the stepping motor, and at the weight used for 
stretching the paper. Any dilatation of the paper strip 
would appear in a difference between the speeds at these 
points. We have used enough small stretching weight to 
avoid any such a difference. 

We monitor the wetting front h{x, t) with a high resolu- 
tion (0.48 million pixels) CCD NTSC camera horizontally 
centered in front of the front plate. We use a personal 
computer to digitize the video image, and to calculate 
the average height h — (h(x,t)) x in real time. We set 
our video system to digitize a 8.4cm wide segment of the 
filter paper strip involving a 110/im (squared) pixel size. 
We control the motion of the paper by a gear mecha- 
nism driven by a stepping motor connected to the same 
computer. According to our calibration, one step of the 
motor corresponds to a paper shift of 5y — 85/im. 

The water front moves upward in the paper due to the 
effective driving force e, which is mainly determined by 
the balance between capillarity and gravity. We maintain 
this driving force constant by holding h at a predefined 
value H. If h exceeds H by Sy, then we pull the paper 
down by Sy. This negative feedback, apart from fluc- 
tuations, prescribes a constant average speed V for the 
paper for a given H . Our measurements at 14 different 
values of H fit remarkably well the power law 



control parameter, H w [2mm, 40mm], which spans ap- 
proximately 2.5 decades of the capillary number. 



V <x H 



-n 



(4) 



where O = 1.594 ± 0.007, and the quoted uncertainty 
is the standard deviation from the fitted value. The re- 
lation (0) holds in the entire investigated range of the 
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FIG. 2. Snapshots of the evolving interface for 
different velocities V. Units are in pixels, (a) V — 
4.03 x 10- 2 cm/sec, (b) V = 1.32 x 10~ 2 cm/sec, (c) 
V = 1.40 x 10 _3 cm/sec. The sampling rate is the same 
constant for all pictures and the artificial shift k between 
the contour lines is proportional to this constant. 

During the experiment, the computer saves the contour 
lines of the digitized images for later analysis. Snapshots 
of the evolving interface at different times are superposed 
in Fig. 2 for three different values of the velocity. We find 
only a few and negligible fl(]|| overhangs, and when they 
occur we remove them in the usual way pjj. 



2 




0.0 1.0 2.0 3.0 4.0 

log (t) 

FIG. 3. The differential autocorrelation function Y 
vs. time is shown by circles. The dashed lines present 
the saturation regime of the autocorrelation function. 
The speeds for the 3 curves are the same as in Fig. 2. 
Larger speed implies larger intrinsic width (intercept), 
and smaller saturation value. 

Using the contour lines we calculated C(t), which has 
less uncertainty than the standard deviation of the inter- 
face height, and is therefore more useful for calculating 
the true exponents |12| . To check the scaling behavior of 
C(t), we form the quantity 

Y 2 (U) - C2( ;g°'g fe) x tf (5) 
h+i h 

where fi,ij+i are successive times. If C(t) oc t@, then 
Y(t)\pr = f3 must scale the same way, but provides a 
more precise (self-consistent) method to determine (3 by 
eliminating any existing additive "correction-to-scaling" 
terms Figure 3 presents our results for the same 

data sets as shown in Fig. 2. We find Y(t) oc 1r for early 
times and this power-law behavior extends over slightly 
more than three decades at the minimum driving force 
(V = 1.4 x 10 _3 cm/s). Within the experimental un- 
certainty, is independent of our control parameter H, 
yielding an average value 

/3 oxp = 0.56 ± 0.03. (6) 

This result can be compared to 2 dimensional theoret- 
ical results and simulations because the small thickness 
of the filter paper (which is in the order of microns) does 
not permit any coarse-grained roughening perpendicu- 
lar to the paper strip. Comparing /3 cxp to ([3a]), we may 
exclude the dominance of annealed randomness, which 
prescribes a much lower value for 0. Moreover, the value 
°f 0qew is also inconsistent with our result. Nevertheless 
quenched disorder with non-linear dynamics produces a 
remarkably close exponent, /3qkpz = 0.6 ||, which sug- 
gests that quenched randomness and non-linearity play 
important roles in the dynamics of vertical imbibition. 



Remarkably, we find no crossover for j3 in the scaling 
regime, as might have been expected in the case of QEW 
IIH, where the temporal scaling is described by Pqew 
for early times and by /3ew for late times (at an inter- 
mediate driving force). As the driving force increases, 
the crossover point t* becomes smaller and the scaling 
region described by /3ew enlarges. From Fig. 3, we ob- 
serve crossover neither as a function of driving force, nor 
as a function of time — even at the at the maximum driv- 
ing force (V — 2.1 x 10 _1 cm/s). The absence of such 
a crossover also indicates the importance of the non- 
linearity in our experimental system. 

It has been observed pljj] that numerical results fall 
into two groups depending on the origin of the non-linear 
term X(Vh) 2 . Kinematics produces a A which vanishes 
at the threshold and the resulting interface belongs to 
the same university class where the non-linear term is 
absent, A ff = 0. On the other side the QKPZ univer- 
sity class is characterized by a non-vanishing non-linear 
term leading to (A e g ^ 0). What is the origin of such a 
non-linearity which leads to /3 cxp « /3qkpz in our exper- 
iment? Recently Tang, Kardar, and Dhar demonstrated 
|]l6f that anisotropic depinning yields a non-zero A c ff at 
the depinning transition. To illustrate the existence of 
such an anisotropy in our system, we consider the filter 
paper as an interconnected network of different capillary 
tubes. The inhomogeneity of this network can be consid- 
ered as a random field with amplitude A. This random 
field is correlated isotropically in space within a distance 
a, but the driving force (pressure drop) is different along 
horizontal and vertical tubes, Fh = / and F v = f — pga 
respectively, because the driving force must compensate 
the weight of the liquid in a vertical tube (p is the den- 
sity of the liquid). Due to this anisotropy, a segment of 
the interface can easier be pinned vertically than horizon- 
tally. Therefore a slope-dependent effective driving force 
F(Vh) — F(Wh) — F c is generated under coarse graining, 
which explains the absence of large overhangs. Tang et 
al. pointed out jl6j that an expansion of F(X7h) around 
its minimum yields to a term A e ff(Vft.) 2 which remains 
finite, independent of V, which is the hallmark of the 
QKPZ universality class. 

Figure 3 demonstrates two effects of increasing the 
driving force, namely the intercepts increase and the 
saturated values decrease. In order to incorporate this 
behavior into the scaling formalism, we assume that 
Eq. (|l|) is valid with the rescaled time and space vari- 
ables t' P = t^V 9i and L' a = L a V~ 8L . As a consequence 
of this assumption, C(t)i,,v scales as 

C(£)l,v ~ V-° L L a T ^L-«// 3 y (et+ ^ )//3 ) . (7) 

This extended dynamic scaling (EDS) hypothesis is fully 
compatible with Eq. (jlj) for a given V, but also describes 
the characteristic role of the driving force in presence of 
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quenched disorder. In the limit of small driving force 
and early times (t [3 V 0t+BL < L a ) the EDS hypothesis 
reduces to C ~ t^V 6t . For t ^> t c , where t c is a critical 
time defined by t c = L a / f3 V-( 9t+9L ' )/l3 , the autocorrela- 
tion function is independent of time and saturates at a 
constant value C ~ L a V L . This asymptotic behav- 
ior is consistent with Fig. 3, as it is expected. In or- 
der to perform a quantitative check of Eq. (Q), we plot 
Y' = CL- a V 6L vs. X' ee t L- a /Pv {et+e ^/ fi for different 
values of Ot and Or,- The best data collapse is found with 
t = 0.37 and 9l = 0.48. Using these exponents we form 
the scaling plot for measurements at 5 different values 
of H and times spanning three decades. From the data 
collapse in Fig. 4, we conclude that our measurements 
are fully consistent our scaling Ansatz. 
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FIG. 4. Scaling plot for five different data sets. 
In addition to the speeds of sets on Fig. 2 two addi- 
tional sets are also drawn for V — 2.50 x 10 _3 cm/sec) 
and V = 4.20 x 10~ 3 cm/sec; Y' and X' are defined 
in the text. The best result is achieved with exponents 
13 = 0.56, 6 t = 0.37, L = 0.48. According to the scal- 
ing function T the different data sets collapse on two 
straight lines, with slopes (3 = 0.56 and zero for early 
and late times respectively. 



our system. The relation h ~ V~ x / n points in this direc- 
tion, although it is clear that Eq. (0) must break down at 
the pinning transition, where V = 0. If we now replace 
£ by V~ v l e in w(e ^ 0,L,t), and compare the result to 
Eq. (0), we see that there is a mapping between the two 
relations with a' = a+OOhjv and z' = z+0(0 t +0L)/((3v). 
Kertesz and Wolf demonstrated that a' and z' are the 
characteristic exponents at the transition point between 
two different morphological phases. It remains an inter- 
esting open question if these exponents have the same 
significance at the pinning transition too, and this ques- 
tion requires more experimental investigation. 

To summarize, we have studied the temporal behav- 
ior of roughening interfaces during vertical imbibition in 
porous media. The height-height autocorrelation func- 
tion C(t)r,y of the meniscus exhibits temporal scaling 
without crossover. The corresponding exponent is inde- 
pendent of driving force, (3 = 0.56 ± 0.03. We conclude 
that nonlinearity plays important role in the dynamics 
of our system. For a compact description of our mea- 
surements, we suggest an extended dynamic scaling hy- 
pothesis, which describes the characteristic role of the 
driving force in the presence of quenched disorder, and is 
fully compatible with the dynamic scaling developed for 
systems with annealed randomness. From the resulting 
data collapse, we find t = 0.37, #l = 0.48. Finally, we 
note that Rubio et al. jl9| measured the interface width 
cr(l,t) versus / in a different experiment, and they re- 
ported a velocity dependent scaling w(l) ~ l a v~ 0A7 , a 
result consistent with our value of Or,. Although our re- 
sult suggests that (3 is independent of the dynamics of the 
system, it remains to be seen whether the moments of the 
quenched noise play an important role in determining the 
universality of exponents. 



In fact there is some theoretical base for the EDS 
hypothesis: Kertesz and Wolf jl7) pointed out that a 
new diverging length £ ~ in the dynamics modi- 

fies the scaling of the interface and £ must appear in the 
scaling form w(e,L,t) ~ £ Q H(t',L'), where V = L/£, 
t' ee t/^ z , and z' ee a' I f3' . Far from the transi- 
tion point H(t',L') - L' a T(t! '/L W/3 ), which involves 
w(e ^ 0,L,t) - i a '- a L a T(tlL z i z - z '). This expres- 
sion resembles Eq. (^), and the similarity is more than 
merely formal. Movement of the meniscus consists of pin- 
ning and depinning processes Jig] . Computer simulations 
demonstrate M that as the driving force e decreases, the 
surface slows down according to V ~ e e and the linear 
size of the pinned regimes £ diverges £ ~ e~ v . These two 
relations can be integrated, with £ ~ V~ v l e . According 
to this picture, we also expect a new diverging length in 
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